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Abstract 

The perturbative renormalization of the Ginzburg-Landau model is reconsidered based 
on the Feynman diagram technique. We derive renormalization group (RG) flow equations, 
exactly calculating all vertices appearing in the perturbative renormalization of the (p'^ model 
up to the order of the e-expansion. In this case, the (^^, ip'^, (p^, and ip^ vertices appear. All 
these are relevant in the sense that they do not vanish at the fixed point. We have tested the 
expected basic properties of the RG flow, such as the semigroup property. Under repeated 
RG transformation Rs, appropriately represented RG flow on the critical surface converges to 
certain s-independent flxcd point. The Fourier-transformed two-point correlation function 
G(k) has been considered. Although the e-expansion of X(k) = 1 /G(k) is well deflned on the 
critical surface, we have revealed a fundamental inconsistency of the perturbative method 
with the exact rescaling of -'^"(k), represented as an expansion in powers of fc at fc — >■ 0. 
We discuss also certain paradox, related to the e-expansion of X(k) above the critical point. 
Apart from the e-expansion, we have tested also a modified approach, where the coupling 
constant u is the expansion parameter at a fixed spatial dimensionality d. Our tests point 
to some internal inconsistency of such a method. 
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1 Introduction 



The perturbative renormalization of the Ginzburg-Landau (or ip^) model has a long history 
(see [Tl[2l[3lH] and references therein). However, as mentioned in [2j, a complete formulation of 
the renormalization group (RG) beyond the order of the e-expansion (where e = A — d>0, d 
being the spatial dimensionality) met mathematical difficulties, which could not be overcome. In 
fact, the e-expansion of the critical exponents beyond the lowest orders is based on an alternative 
approach, which relies on the Callan-Symanzik equation O [U [7] . The latter one represents a 
scaling property of the (^^ model [5] , and the method is based on a set of assumptions [7] . Since 
the perturbative RG theory is not rigorous, stringent tests of its validity and consistency make 
sense. Our aim is to perform such tests. 

We have tested the expected basic properties of the RG flow. Following the idea in [8], we 
have checked the semigroup property i?siS2/^ = -Rs2-^siM' where Rg is the RG operator with 
scale factor s > 1 acting on the set of Hamiltonian parameters fi. We have tested the expected 
s-independence of the fixed point /x*, as well as the scaling of the Fourier-transformed two-point 
correlation function G'(k) on the critical surface. 

Another strategy of verification has been used in [9] , considering a four-dimensional (d = 4) 
model with Gaussian measure modified in such a way to simulate d = 4 — e dimensions with e 
small and positive. This method allows to control rigorously the remainder of the perturbation 
series. The obtained results [9] confirm the existence of the non-Gaussian fixed point at a dis- 
tance 0{e) away from the Gaussian one as proposed earlier [1] by the e-expansion. However, 
the expected decay of the two-point correlation function appears to be canonical (i. e., Gaus- 
sian, see Introduction part in [9j) in disagreement with that provided by the e-expansion and 
observed in ordinary spin systems of dimensionality d < 4 like, e. g., three-dimensional and two- 
dimensional Ising models. In view of these observations, the classical (used in the e-expansion) 
way of introduction of noninteger spatial dimensionality d via an analytic continuation from 
d-dimensional hypercubes appears to be more meaningful than that in [9]. In fact, there is no 
unique definition of the noninteger d. It can be introduced in a more physical way as a suitable 
"fractal" dimension of an irregular lattice \10\ [TT] . 



2 Diagrammatic formulation of the renormalization 

As a starting point, we consider a 99^ model with the Hamiltonian H defined by 

H/T = j"^ ^\^) + |(V(p(x))2 + 111 ¥^'(xi)n(xi - X2)v^\x2)d^idK2 , (1) 

where the order parameter 9?(x) is an n-component vector with components v'i(x), depending 
on the coordinate x, and T is the temperature. The field (pj{x) is given in Fourier represen- 
tation by (/3j(x) = V'^^"^ J2k<AfjM^^^^^ where V = L'^ is the volume of the system, d is the 
spatial dimensionality, and A is the upper cutoff of the wave vectors. The Fourier-transformed 
Hamiltonian obeys the equation 

-f = -^E(- + -k2)l^.kP (2) 

i,k 

- g^"^ X! ¥'i,ki9'j,k2'"ki+k2¥'j,k3 9'i-ki-k2-k3 • 
j,j,ki,k2,k3 
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In the Feynman diagram technique O [12], the second term in ([2]) is represented by a fourth 

order vertex ^ ^ (where field components with vanishing total wave vector are related 

to the solid lines and the remaining factor — to the dashed line). In the following we will 
consider only a particular case ii(x) = u6{x) or Uk = u. The fourth order vertex then can be 
depicted as by shrinking the dashed line to a point (node). We will mostly use this 

simplified notation. One has to remember, however, that two of the vertex lines are related to 
i-th component (i = 1, . . . ,n), and the other two lines — to j-th component of a field vector 

(including the possibility i = j). It is shown explicitly in the representation ^ <^ , where 

these pairs of lines are separated. 

In the exact Wilson's RG equation the scale transformation, i. e., the Kadanoff's transforma- 
tion, is performed by integrating over the Fourier modes with wave vectors obeying A/s < /c < A. 
It is the first step of the full RG transformation. At this step the transformed Hamiltonian H' 
is found from the equation 

i, k/s<k<k 

where ^4 is a constant. Note that (/9j^k is a complex number and (/7i,-k = V'i k holds (since '/'i(x) is 
always real), so that the integration over (/9j^k in © means in fact the integration over real and 
imaginary parts of (/?j^k for each pair of conjugated wave vectors k and — k. In practice Eq. ([3]) 
cannot be solved exactly. It is done perturbatively, as described in [2]. The perturbation terms 
can be found more easily by means of the Feynman diagrams. 

In the perturbative approach, Hamiltonian is split in two parts H = Hq + Hi, where Hq 
is the Gaussian part and Hi is the rest part considered as a small perturbation. The first 
and the second term on the right hand side of ([2]) can be identified with —Hq/T and —Hi/T, 
respectively. Since the term with r also is considered as a small perturbation (r ~ e holds in 
the e-expansion in d = 4 — e dimensions) , it is suitable to include it in Hi . In this case diagram 

expansions are represented by vertices ^ ^ and — • — , where the latter second-order 

vertex corresponds to the term with r. 

It is convenient to normalize Eq. ^ hy Zg = J exp(— ifg/r) Hi, A/s<fe<A '^V'i.kj where Hq is 
the part of Hq including only the terms with A/s < k < A. We have In Zg = As L'^ at L ^ oo, 
where Ag is independent of L, and the normalization yields 

- (H'/T) + iA- Ag)L'' = -{H'JT) + ln(exp(-/?i/T))o , (4) 

where (•)o means the Gaussian average over the field components with A/s < < A, whereas Hq 
is the part of Hq including only the components with k < A/s. Like free energy, ln(exp(— //i/T))o 
is represented perturbatively by the sum over all connected Feynman diagrams made of the ver- 
tices of —Hi/T by coupling those solid lines, which are associated with wave vectors obeying 
A/s < k < A, according to the Wick's theorem. A diagram can contain no coupled lines. There 
is also a set of diagrams with all lines coupled. The latter diagrams give a constant (independent 
of 93i,k) contribution which compensates the term (A — As)L'^ in dH). The other contributions 
correspond to —H'/T. 

Each term of ln(exp(— //i/T))o comes from a diagram (or diagrams) of certain topology 
and is given by a sum over wave vectors which fulfil certain constraints such that k < A/s 
holds for uncoupled (external or outer) solid lines, associated with field components (pi^k, and 
A/s < k < A holds for coupled solid lines. Besides, the sum of all wave vectors coming into any 
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of the nodes is zero, and the same index i is related to the sohd hnes attached to one node of 

any vertex ^> ^ or — • — . The Gaussian average Go(k) = {^Pi,k.^Pi~'k)o = (I ¥'j,k |^)o 

is related to a coupling line with wave vector k. Here Go(k) is the Fourier transform of the 
two-point correlation function in the Gaussian approximation. If only the term with k^ in ([2]) 
is included in Hq/T, then we have Go(k) = l/(c/c^). Including also the term with r, we have 
Go(k) = l/(r + cA;2), but in this case we need finally to expand Go(k) in terms of r. The second 
method yields the same results as the first one: such an expansion generates the same terms, 
which are obtained in the first method by extending coupling lines in the diagrams originally 
constructed without vertices — • — to include all possible linear chains made of these vertices. 
We will use the first, i. e., more diagrammatic method. Note that, when the renormalization 
procedure is repeated, a continuum of additional vertices appear in the expansion of —Hi/T, 
including ones with explicit wave-vector dependent factors related to the solid lines, which then 
also have to be taken into account. 

In fact, the first step of RG transformation implies the summation over wave vectors of the 
coupled lines in the diagrams, resulting in a Hamiltonian which depends on c^j^k with k < A/s. 
The RG transformation has the second step [2] : changing of variables k = sk and rescaling the 
field components (pi^]^ — )• where rj is the critical exponent describing the ~ k^'^^'^ 

singularity of the Fourier-transformed critical two-point correlation function at A; — >• 0. The 
upper cutoff for the new wave vectors k is the original one A, whereas the density of points 
in the k-space corresponds to ■s'^ times decreased volume. Therefore we make a substitution 
V = s~'^ V. In the thermodynamic limit we can replace by ^ consistently increasing the 
density of points in the wave vector space. Finally, we set k — t- k and obtain a Hamiltonian in 
original notations. 



3 Renormalization up to the order of £^ 
3.1 RG flow equations 

Here we consider RG fiow equations including all terms up to the order O (e^) , starting with the 
renormalization of the coupling constant u. In fact, this is the lowest order of the theory, since 
it allows to find the fixed point value u* up to the order of e. Our aim is to show that we can 
easily recover the known results using the diagrammatic approach described in Sec. [2j It allows 
also to write down unambiguously (i. e., without intermediate approximations) the formulae for 
all perturbation terms. 

Important statements here are that the renormalized values of r and u are quantities of order 
0{e). Besides, rj = O [e"^) holds within the e-expansion and, for any finite renormalization scale 
s, the variation of c is of order O (e^) [2]. Hence, performing the RG transformation Rs with 
a finite s, the only diagrams of ln(exp(— ffi/T))o contributing to the renormalized coupling 
constant up to the order of are and /<C^^ (this statement can be verified in 

detail based on a complete renormalization discussed further on). The first one provides the 
original ip'^ term in Eq. ([2]), which is merely renormalized by factor s^~'^'^ at the second step 
of the RG transformation. The second diagram, which is constructed of two vertices , 
yields 

>0< S^"^'' V'^ '^i,ki¥'i,k2 9'i,k3V'i-ki-k2-k3 X 

i,i,ki,k2,k3 

X [(4n + 16) Q((ki + k2)/s, s) + 16 Q((ki + ks)/s, s)] . (5) 
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Here Q is given by 



Q(k,s) = -1^ / r'|k-qr2^(|k-q|,s)A, (6) 



A/s<(j<A 

where J-{k,s) = 1 if A/s < < A, and J^{k,s) = otherwise. To obtain this result, we have 
deciphered the ^<Z>\ diagram as a sum of three diagrams of different topologies made 

of vertices ^> <^ , i- e., ^ — O — ^ ) ^ C ! , and | j , providing the 

same topological picture ^<Z>\ when shrinking the dashed lines to points. Note that any 

loop made of solid lines of ^ <^ gives a factor n, and one needs also to compute the 

combinatorial factors (see, e. g., [21 [12]). For the above three diagrams, the resulting factors are 
4n, 16, and 16, which enter the prefactors of Q in ([5]). 
Quantity Q(\<./s,s) has a constant contribution 

Q{0,s) = KdA-'{s' -l)/e , (7) 

as well as a k-dependent correction A{k,s) = (5(k/s,s) — Q{0,s) vanishing at A; = 0. Here 
Kd = S{d)/{2^Y, where S{d) = 2t:'^I^ /T{d/2) is the area of unit sphere in d dimensions. 

The term provided by the Q(0,s) part of Q(k/s,s) in (|5]) is identified with the (p^ vertex, 
contributing to the renormalized coupling constant u' . It is consistent with the general form 

XI Q(ki,k2,k3,s,n,e)93i,ki<^j,k2<^j,k3V'i-ki-k2-k3 , (8) 

ij,ki,k2,k3 

of quartic {^^) terms, where the contribution corresponding to the ordinary (p^ vertex is uniquely 
identified with one provided by the constant part Q (0, 0, 0, s, n, e) of the weight function Q. 

Taking into account also the contribution coming directly from vertex, the RG flow 

equation reads 



2 Kd{n + 8) 
U - U X 



2c2A^ 



+ 



O (e') , (9) 

where u' is the renormalized and u is the original coupling constant. The expansion in e yields 
the well known equation [2] 

u' = u + euhis - V? B his + O (e^) , (10) 

where B = K^^n + 8)/(2c^), with the known fixed-point value at u = u* = B^^e + O (e^). 
The correction term A(fe, s) changes the renormalized Hamiltonian as follows: 



2 

, J, .v^ X 

i,j,ki,k2,k3 

x[(l + n/4)A(| ki+k2 \,s) + A{\ ki+ks \,s)] . (11) 



It can be represented as a sum of two fourth-order vertices including prefactors A(| ki + k2 |, s) 
and A(| ki + ka I, s), respectively, related to two of the vertex lines. 

In the following we will introduce a more suitable diagrammatic notation. In the one- 
component case n = 1, the term appearing in ([5]) is represented as 

.2-2^ >::::::< (12) 

16 
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where is the Feynman diagram in which the Gaussian propagator l/(cA;^) with 

k £ [A, sA] is related to the dotted couphng hnes. In other words, the propagator is multiplied 
with the cut function J^(k/s,s). In general (also for other diagrams considered in our paper), 
the field components supplied with the factor y^~^<- are related to 2m external lines, the rest of 
such factors being absorbed in the k-space integrals. An additional factor in ()12p comes from 
the rescaling of wave vectors in such a way that the integration now takes place over /c G [A, As] 
for the internal lines. This procedure is necessary in order to represent the result as a Feynman 
diagram with vanishing sum of the wave vectors entering each node. 

Subtracting from the k-space integral its value calculated at k = for the external lines, 
the result is represented as 

lb 

Similar notation is applied throughout the paper, where the line crossing the diagram indi- 
cates the subtraction of zero~k contribution. The diagram representation of this kind has been 
proposed in [13]. In the general n-component case, we denote by 

. ( >:^< ) ^ i >-^-< . i >-^/L - I (") 

the sum of all such kind of diagrams, which reduce to when the dashed lines shrink 

to points. In (|14p the coefficients are normalized in such a way that their sum is 1 at n = 1. 
Note that in the diagram technique with such vertices, one has to count closed loops made of 
all kind of lines related to the propagator l/(cfc^). In this notation (llip becomes 

H/T ^ (H/T) - s^'-^^u^ >■■■/:■< ) (15) 

Analogous diagrammatic representation is possible for all vertices appearing in a repeated renor- 
malization. In accordance with such notation, the ordinary vertex in ([2]) will be represented 

as -{u/8) > < . 

For a complete renormalization up to the order of e^, one has to take into account also other 
contributions of this order, including the vertex 

s(> <) ^ >-J L_< 

= C"V"^ V'i,ki¥'i,k29'i,k3'/'j,k4V'«,k5'/'i,-ki-k2-k3-k4-k5 (16) 

ij,Z,ki,k2,k3,k4,k5 

X I ki + k2 + kg T{\ ki + k2 + kg I /s, s) 

appearing due to the coupling of two vertices , as well as the quadratic term of the general 

form (1/2) ^0(k) | ipi^^ P provided by diagrams with two external lines. In the representation 

i,k 

0(k) = r + ck^ + 0(k), it corresponds to the ip'^ vertex in ([2]) with one additional term. 

In the rest part of this section we will consider the structure of the renormalized Hamiltonian 
and the related RG flow equations, providing the proofs afterwards in Sec. 13. 2i Thus, the 
renormalized Hamiltonian has the form 

f = ^E(- + ^k2 + ^(k))l¥'.kP+^>--< 

i,k 

+ a4 S ( >■/'.'■< ) + ae S ( > < ) + O [e^) . (17) 
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Here the wave vectors of the dotted Hnes are in the interval [A, (^A], where the parameter C has 
the initial value 1 and is transformed according to 



C' = sC 



(18) 



in a process of repeated renormalization. Here C is the previous and C is the new (renormalized) 
value of this parameter after the current RG transformation (RGT) with the scale factor s. The 
weight coefficients 04 and uq in this order of the e-expansion obey simple RG flow equations 



9u^ 

„,2 



+ 



+ e 



(19) 
(20) 



relating the new values of these parameters in each RGT to the previous value of u. The initial 
values of 04 and og are not important, since the corresponding diagrams vanish at = 1. The 
RG flow equations for the Hamiltonian parameters u and r read 



n + 



u 



-u 



s 

n + 2 



—u oU 



^ n + 2 



u qO 







16 



04 o(/) -24a6 0© 



+ 0{e 



(21) 



(22) 



where o"^^^ is the diagram with amputated four external lines having zero wave vectors, and 
the other diagrams are defined analogously. The wave vectors of the internal solid lines are in 
the range [A/s,A], whereas those of the dotted lines — in the range [A, A(^]. Eq. (I2ip is the 
same as Q, only the irrelevant for this order of the e-expansion factor s"^'' is omitted, like 
also the factor in (j22p . Taking into account (|19p - (j20p and the fact that u is renormalized 
only by O (e^) in one RGT, the diagram expression in the last line of ()22p can be written as 

-\{n + 2)v? 







+ 3 o(;/:) + 3 



The parameter c enters all diagrams. In the actual order of the e-expansion its variation does 
not produce additional terms in (f2T]) and (p2]) . However, one has to update its value according 
to the equation 



n + 2 



n + 



where 



with 



H 



lim 

fc-*-o 



D{k) - D{0) 



D{k) 



+ 3 kt— ) + 3 k' 



(23) 



(24) 



(25) 



Here the symbol k indicates the wave vector related to the amputated external line. Finally, 
the RG flow equation for the function ^(k) reads 



9'{k) 



e(k/s)-^^ik/s))+0{e') 



(26) 
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where 

$(k) = Dik) - D{0) - k'^U . (27) 
Using the diagram technique, the function 0(k) can be represented as 

^(k) = -^n^ k(|) +0(e^) , (28) 

where the wavy hne indicates that both the constant contribution and that proportional to fc^ 
(with the proportionahty coefficient determined at /c — )• 0) are subtracted. Besides, the range of 
the wave vectors for the internal lines is [A, A^], where C, is determined after the actual RGT. 

3.2 Proof of the RG flow equations 

In Sec. 13. li the RG flow equations are given without proof. Their proof consists of a verification 
that, at each RGT, the Hamiltonian keeps the form ()17p in accordance with the given update 
rules for all parameters. Consider first the 93^ vertex at n = 1, summing up all diagrams of this 
topology, which appear in H/T at a given RGT: 

= -y > < +0(e^)=a',^ < +0{e^) . (29) 

Here the wave vectors within [A, As], [As, AC'], and [A, A^'] correspond to the solid, dashed, 
and dotted coupling lines, respectively, where C' = sC,. (These dashed lines should not be 
confused with those of ^ ^ .) In ([29|) . the vertex ^ — ^ is produced by coupling 

(and rescaling) two vertices ^^>\ , whereas ^ ^ comes from the rescaling of already 

existing before the given RGT cp^ vertex. We have omitted the irrelevant for this order of the 
e-expansion rescaling factor resulting from the renormalization of parameter c. We have used 
the substitution oe = — n^/8 + O (e^), which holds according to the assumption that our RG 
flow equation (f20]l was satisfied and u was changed only by O (e^) in the previous RGT. It is 
true for the first RGT, so that the relations (j29]) provide the proof of (f20|) by induction. It is 
obviously valid for any finite number of RGT. 

The generalisation to the n-component case is trivial here, since there is only one way how 
to couple two vertices ^ ^ to form the sixth-order vertex. 

Similarly, the summation of terms provided by the diagrams of ^<Z>\ topology, sub- 
tracting the contribution included in the ordinary vertex, gives us 

s'^-'^ (-^ >0< +a, >7t< +9a6 >ZZK ) 

= -^(>;zK + >'/-< +2>szK)+o(.^) 

= >7^< + O {e') = a', >:t< + O (e^) (30) 

with the same k intervals for different coupling lines as in (j29p . It proves (jl9p for a finite 
number m of RG transformations at n = 1. Here the three diagrams in the first line of (j30p 
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come from two vertices, rescaling of the already existing vertex with the coefficient 04, 

and from the 99® vertex via coupling its two lines. Besides, holds, since 

the diagram ^<C^>^ vanishes if the external lines have zero wave vectors. This derivation 
is easily generalised to the n-component case represented by (jl7p : the same relations, only 
multiplied with the corresponding weight factors in (I14p . hold for diagrams of each topology, 

when is deciphered as the diagrams on the right hand side of (|14p . 

A subtle question here is whether the small corrections, omitted in (I29p and (I30p at each 
RGT, do not lead to larger than O (e^) deviations from the assumed form of the renormalized 
Hamiltonian when the number of RG transformations m growth unlimitedly. A rigorous answer 
to this question can be found, comparing the results of the actual update rules with those 
provided by a renormalization scheme in which all these corrections are taken into account. 
Such a scheme will be considered in Sec. [H Here we note only that the actual RG flow equations 
are valid within the given accuracy for any m. 

The derivation of ([2T]) has been already discussed in Sec. 13. li Recall that the updated value 
of the coupling constant u of the ordinary tp^ vertex is related to the constant part of the weight 
function in ([8]), which comes from the rescaled already existing ordinary ip^ vertex, as well as 
from other diagrams with four external lines generated in the actual RGT. The coupling of two 

lines of the vertex ^ ^ gives vanishing contribution here, so that the coupling ^^>C!^>\ 

is the only relevant one at the order of e^. It results in (j2ip . 

Similarly, the updated value of the parameter r is identified with the constant part of the 
weight function 0(k) in the general representation of the quadratic term (1/2) ^0(k) | y^j k 

and is generated by the diagrams with two external lines. It results in (|22p . where the terms with 

04 and ag appear via coupling the lines of the vertices and ^ ^ , respectively, 

whereas other diagrams in (j22p are produced by the vertices and — • — . 

The RG flow equation (j23p for the parameter c is similar to that for r with the only difference 
that here we single out the contribution to ^(k), which is proportional to k^. Subtracting both the 
constant contribution and that of k"^, we arrive at the equation for 0(k.) (126p . Its diagrammatic 
form ()28p is proven by summation of diagrams like, e. g., in Eq. (j30p . 

Finally, we check that all terms up to the order of are already included. In particular, the 
part of ()17p with factor ^(k)/2 can be seen as the vertex — • — supplied with such a weight 
factor. However, its coupling to other vertices gives only a contribution of order O (e^). Recall 
that we do not consider the constant part of Hamiltonian, which is independent of the field 
configuration and is represented by closed diagrams without external lines. 



3.3 Estimation of k-space integrals 

It is important to know whether the k-space integrals over the internal lines of the diagrams, 
appearing in our representation of the renormalized Hamiltonian, are convergent when the RG 
transformation is repeated unlimitedly many times, i. e., at C — ^ 00. 
First, we note that the diagram o'------* iii four dimensions, i. e., 

1 f d^q K4 fdq 

A<<j<AC A 

diverges when the region q G [A, A^] of wave vectors for the internal lines is extended to [A, 00] 
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at C — ^ oo. Similar diagram, where the zero-k contribution is subtracted, 

i{kx) = k--t-- = k-::::> - o-::::> (32) 

with k corresponding to the sum of wave vectors for the two amputated external lines, is con- 
vergent at d = 4, as well as in d = 4 — e dimensions, where we have 

A<fc<AC \ ' ' / 

} ( P i^q^ + 2kqcose + k'^,C] i\ 

c^hiner-^de I {[ l' + ^kqcose + k- k-)' 



with J-'{k, (^) = 1 if A < A; < A(^, and J-{k, Q = otherwise. The last identity in (j33p is true for 
any integer d, according to the well known formula for integration in spherical coordinates, so 
that we can put here e as a continuous parameter in the usual sense of the e-expansion. 

In the following calculations of this subsection, we set c = 1 and A = 1 for simplicity. 
Considering the large-g asymptotic at ^ = oo, we can set = 1 in (|33p . which gives correct 
result up to a term vanishing at q ^ oo. Further on, we use the decomposition 

^ ^=/o(A:,(?)+A/(k,q), (34) 



fc2|q + k|2 fc4 
where 

An important property of A/(k, q), defined by ([M]) and ([35]) . is 

J A/(k,q)d^A; = at q ^ oo . (36) 
fc>i 

It is true because 

A:>1 k>l 

holds at g — )• oo. The —lnq + 1/2 asymptotic is evident for the second integral. For the first 
one, it is verified by integration over k (via the substitution p = k + q cos 6) and then - over the 
angle 9 in spherical coordinates. According to the decomposition (I34p . we have 

/(g,oo) = j^rfj^- ' I +A/(g) at q ^ oo , (38) 

where the term Kd{-} is the contribution of fo{k,q), whereas A/(g) comes from A/(k, q) and 
is equal to 

OO TT 

A/(g) = ;^r[(3-£)/2] I ^''"'"^^ / ^fWq,0)(.sm9f-'de , (39) 
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where 



Af{x,0) 



1 + (1-2;^) 0(1 -x) . 



1 + 2x cos 9 + 

Here G(x) denotes the theta step function, and the known identity 



r(^) 



(40) 



(41) 



has been used to obtain ([39ll . According to (1361) . AI{q) vanishes at d = 4. Using this property, 
from (1391) we obtain 



2 

AI{q) = -Kd A (-£ + In q) +c£^ + (e^) , 



where 



Inx 



00 TT 

dx 



X 





A= / dx Af{x,e)sm^ede + / — / A/(a;,0)sin^e ln(sin^)d^ 



1 2A 
4 vr 



(42) 



(43) 



and c is a constant at g — )• 00. According to (j38]l and (142]1 . we have 
/(g, 00) = -fCd 

+ — (-21n3g + 31n2gj +ce 



In g + ^ + I ( In^ g - In 5 ) + 



+ e 



e — e In g 
at (7 — >• cx) , 



(44) 



where c is a constant. 

In a more general case, for /(g, C,) we obtain 

/(g,C) = %,oo) + W(g,C), 

where 51{q, Q) has the asymptotic form 

W(g, C) = q^'fiC/q, e) at g, C ^ 00 

where 



(45) 



(46) 



f{y,e) 



Kdm-e)/2] 
0Fr[(3-e)/2] 



X ^ ^dx 



x'^ e(l + 2x cos 6* + - 
1 + 2x cos 6 + x"^ 

-1) (sin0)2-^de 



X ^ ^dx 



1 + 2x cos + x^ 



(sine)2'"d0 



(47) 



Based on these results, we can evaluate the large-q (i. e., C < g < C for large enough C at 
a given k) contribution to the diagram kv/ '-^ iii four dimensions as 



(2vr) 



i-^ 1 [/(|k + q|,C)-/(g,C)]9"'A. 



(48) 



C<g<C 
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It, in fact, diverges at — )• oo, and the divergent part comes from the — Inq term in I{q, 00), 
providing the contribution 

Kl f , f , f 2kqcose + o ^ ,^ 
qdq / In 1 + — 2 sin^ 6 dO . (49) 

^ c ^ ^ ^ 

Using the expansion in powers of k/q, we see that (|49|) logarithmicahy diverges at ^ ~^ 00. 
However, if we subtract the term oc k'^, the result is convergent. It implies that the diagram 

appearing in has finite value at C — ^ 00. 

.....q 

Similarly, we can treat the large-wave-vector contribution to the diagram k><i^ further 

considered in the renormalization up to the order of e^. This graph contains the subtraction of 
the zero-vector contribution from the internal block , as well as from the whole diagram. 

It is defined (at c = A = 1) by 

=J^ J {c'|qi+kr'-F(|qi + k|,C)/(|q-qi|,C)-gr'%l,C)}Al 

1<91<C 

(50) 

in four dimensions. Using the asymptotic estimates of I{q, 00) and SI{q, (), we find that the tail 
of large qi is convergent for any given k and q. 

The actual integrals (I49p and (I50p have been estimated in four dimensions (d = 4). Those 
ones, which are convergent at d = 4, are convergent also in d = 4 — e dimensions for small 
(positive and also negative) e when considering the k-space integrals as continuous functions 
of d in the usual sense of the e-expansion. It is because the contribution of large wave vectors 
changes only slightly. 



3.4 Testing the semigroup property 

Here we test the semigroup property 



(51) 



discussed already in Sec. [TJ In fact, we verify that the renormahzed Hamiltonian after two 
subsequent RGT with scale factors si and S2 is the same as after one RGT with the scale factor 

siS2- It can be seen most easily for the vertices ^ ^ and , as consistent with 

Eqs. (jlSp to (j20p and the fact that u is renormahzed only by an amount of O (e^) . The semigroup 
property for the coupling constant u easily follows from (I2ip . After the first RGT we have 



n + 8 2 
u — u 0' 



+ 0[e 



(52) 



where the dotted line refers to the k interval [A/si, A]. After the second RGT the value of the 
coupling constant, denoted as m", becomes 



si 



u — u s^^ O"'--" 



+ 



(53) 
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where the dashed hne refers to the interval [A/(siS2), A/si]. Here the rescahng of the wave 
vectors from k G [A/s2,A] to k G [A/{siS2), A/si] has been performed for convenience. Insert- 
ing (j52|) in ([53l) . and taking into account that u = O (e), we obtain 



u 



{siS2f U- 



n + 8 
2 



u^l Q<.'.'.'-' + 



{S1S2Y u 



n + 8 
2 



(54) 



where A; G [A/(siS2)) A] corresponds to the sohd couphng lines. It is obvious that (|54p is identical 
to the result of one RGT obtained from (|21|) with s = S1S2, which proves the semigroup property. 
Analogous diagram summation proves this property also in the next order of the e-expansion |13] . 

The semigroup property for the parameter r, when performing the RG transformations of 
the initial Hamiltonian with Q = 1, is proven straightforwardly by the same method. Namely, 
the result of the second RGT (applying ([22]) ) is represented by the diagrams containing the 
dotted and the dashed lines with k S [A/si,A] and k G [A/(siS2), A/si], respectively. These 
diagrams sum up to give ones with k E [A/ (S1S2), A], which correspond to those of the one-step 
renormalization with s = siS2- In distinction to the case of u, here more complicated terms 
appear, including different lines in one diagram. The situation is less trivial when starting 
with Q > 1. In this case all diagrams can be decomposed in such ones, which contain three 
types of lines with k G [A,AC], k G [A/si,A] and k E [A/(siS2), A/si], respectively. Then we 
prove the semigroup property by checking that such diagram representation for the two-step 
renormalization is identical (up to the considered order of the e-expansion) to that for the one- 
step renormalization. The semigroup property for the parameters c and 0(k) is proven by this 
method, as well. 

3.5 The fixed point 

The RG flow described by the truncated perturbative equations of Sec. 13.11 has certain fixed 
point as a steady state solution of the RG flow equations. We shall mark the fixed-point values 
by an asterisk, except the value of c, since the fixed point exists for any given c, as it will be 
seen from the following analysis. The fixed-point value of the coupling constant u, i. e., 



has been already mentioned in Sec. 13.11 According to ()19p and (I20p . we have al = -(9'u*^)/16 + 
O (e^) and Og = — u*^/8 + O (e^). It is generally expected that the fixed point is independent 
of the scale parameter s. In the exact renormalization, it is a consequence of the semigroup 
property. In the perturbation theory, it is expected to hold at each order of the expansion. 
The above considered fixed-point values obey this requirement within the given accuracy. It is 
obviously true also for the r* value in the lowest order of the expansion, i. e., 



of the fixed-point Hamiltonian. Therefore, we have performed some tests. Using (155p . (I19p . 




(55) 




(56) 



However, it is less obvious for the next-order correction to r*, as well as for other parameters 
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and (j20p . we find from (\22h that r* is independent of s within the error of O (e^) if u* is 
s-independent up to the order and 

/(s) = --^^ In s + const -(s^-l) (57) 

holds, where 



fis) = s'( oG +3o(/) +3 00) • (58) 



Here the solid lines have wave vectors within k £ [A/s,A], whereas the dotted lines - within 
k G [A, A(^] at C — ^ CO, and all terms are evaluated at d = 4. On the other hand, we prove that 
the function ([58]) obeys the equation 

slf{si) + f{s2) = f{siS2)-3oO X 0-:::::> x{siS2)\ (59) 

where the dashed lines refer to the interval k G [A/(siS2), A/si] and the dotted lines — to the 
interval k G [A/si,A]. The proof consists of a straightforward checking of the corresponding 
diagram identity by the decomposition described at the end of Sec. 13.41 Note only that the wave 
vectors in the diagrams of /(S2) are rescaled by the factor As a test of consistency, we 

verify that ([FT]) is a solution of ([5^ . 

As regards the parameter c, it appears as a common factor in the fixed-point equation 
obtained by setting c' = c and u = u* in (|23|) . taking into account that u* cx c^, whereas the 
diagrams involved are proportional to l/c^. Hence, if this equation is satisfied, then it holds for 
any c. In fact, the equation c' = c reduces to one for the exponent rj: 



2{n + 2) ~ , , ^2 

// ill i; = — - 

where 



is calculated at d = 4 with q £ [A/s, A] corresponding to the solid lines and q G [A, AC,] — to the 
dotted lines in the diagrams. The expected universality of the critical exponent r] follows from 
the diagram identity 

n(si) + n(s2) = n(siS2) , (62) 

which implies that n(s) oc Ins holds and, therefore, r/ is independent of s. The identity (|62p 
is proven by the same method as (j59p . In principle, the critical exponent rj can be calculated 
from ()60p at any given s. However, from a technical point of view, it is convenient to consider 
the limit s — )• oo, since only the first diagram in ()6ip provides the singular contribution ~ Ins 
in this case. Hence, we have 



r? = — - — ^ — ^ „ X lim lim / „ , 



=1 



+ 0(e^) . (63) 



Rescaling the wave vectors by factor s/A, lim lim j-j in ( 1631) reduces to the expansion coefficient 
at A:^ Ins in the asymptotic large-s expansion of the diagram kVy^;.' , in which g G [1, s] 



corresponds to the dotted lines. This coefficient can be calculated from (I49|l at C = •s, and is 
equal to —Kl/A. It yields the well known result (see, e. g., [2j) 
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Closing this section, we note that the independence of the fixed-point function 0*(k) on the 
RGT scale parameter s can be easily seen when using the diagrammatic representation (|28p . 
taking into account that this diagram converges to certain s-independent value at C — ^ oOj as 
shown in Sec. 13.31 

4 A refined diagrammatic renormalization 

Here we consider a refined renormalization scheme, taking into account the correction terms 
omitted in the derivations of Sec. 13.21 In this case the only approximation is that the vertices 
of higher than O (e^) orders are omitted. 

To simplify the notation, consider first the one-component case n = 1, where the symbol 

5 in (|17p can be omitted. The dashed lines, similarly as dotted and solid coupling lines, are 
related to the propagator l/{ck'^) in the notations used in this section. 

After m RG transformations with a constant scale factor s, the terms with and in (I17p 
now look as 

£ (04),/ ( >X< ) , (65) 
j,i=i ^ 

m 

EK). ( > < ) . , (66) 

where the indices j and / are used to distinguish between diagrams with different k intervals 
related to the coupling lines. Namely, the dotted lines in the diagrams refer to fc G [As^~^ , As^], 
whereas the dashed line — to /c G [As'~^, As']. In the RG flow equations of Sec. 13.11 we have 
(0-4) ji = ^4 and (ag)^ = clq, whereas here these weight coefficients are slightly different. 

The generalisation to the n~component case is possible by using the decompositions of the 
kind (|55p and (|66p for the diagram of each topology when deciphered as one composed of the 
vertices ^ instead of . 

Within the O (e^) error of the actual truncation scheme, the weight coefficients (04)^7 and 
(ag)^ after m RG transformations, i. e., {ai)ji{m) and (o6)j(m), are given by the recurrence 
relations 

(67) 

2<j<m, m>2 (68) 

(69) 

i(m-l) : 2<j<m,m>2 (70) 
: 2<j<m,2<l<m,m>2 (71) 

where Rm = c{m)/c{m — 1) is the ratio of c values in two subsequent RG transformations. We 
have denoted by u{m) and c(m) the values of u and c, respectively, after the m-th RGT, the 
initial values being given by n(0) and c(0). Eq. (j67p represents the contribution due to the 
coupling like ^ — ^ and rescaling of two ordinary f'^ vertices, whereas the coefficients (j68p 
are the weight factors of the corresponding terms rescaled repeatedly. Eq. ([69]) corresponds to 
the coupling and rescaling of two ordinary ip'^ vertices like ^«d^^ , whereas (|70p — to the 



8 

{ae)jim) = Rms'^''^'^ {a(i)j-i{m - 1) : 

(a4)ii(m) = -Ai?^s2.-2^n2(m-l) 
Id 

9 

(a4)ij(m) = (04)^1 (m) = -ii^s^^"^'' (ae)^- 
{a4)ji{m) = (a4),_i,/_i(m - 1) 
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coupling of two lines of the vertex ^ ^ with k S [As-'"^, As^~^ related to the dotted line, 

followed by the rescaling of the diagram. Although there are 9 combinatorial possibilities to 
couple these lines, the factor 9/2 in ()70p is twice smaller. It is because the corresponding term 
appears twice in the decomposition ([65]) . i. e., as (04) ij and (a4)ji. Finally, ([7T]) comes from the 
rescaling of the quartic terms ([69]) and ([70|). The actual update rules ([67|) to ([7T]) are consistent 
with the fact that the range of the wave vectors for the internal lines of any diagram is rescaled 
by the factor s in each RGT, taking into account the corresponding rescaling factors discussed 
in previous sections. Here, however, additional rescaling factors Rm and appear, which 
have been neglected before. They originate from the renormalization of parameter c: since the 
diagrams are defined so that factors l/{ck'^) correspond to the coupling lines, these factors need 
to be redefined with new value of c after each RGT. Hence, the weight coefficient for a diagram 
with / such coupling lines gets a factor RI^. 

Eqs. ([67|) to ([TT]) can be applied also to the general n-component case. In this case (04)^7 
are the common factors for the weighted sum of diagrams of three different topologies appearing 
in (jl4p with the corresponding k intervals for the coupling lines. It is obviously true because 
the diagrams of these three topologies always appear with certain relative weight factors, which 
depend only on n and are the same as in (I14p . 

Since u and c are varied only by an amount of O (e^) in one RGT, it easy to verify that 
Eqs. dSZ]) to dZH) provide (a4)j7(m) = (a4)ii(m) + O (e^) and {aQ)j{m) = {ae)i{m) + O (e^) for 
any m and j < M, I < M at any large but finite M. Besides the values of these coefficients are 
consistent within the O (e^) error with the corresponding values of 04 and og provided by (jlDh 

and ([20]) after m RG transformations. Since the diagram is convergent when the 

wave vectors within [A,c«] are related to the dotted lines (see Sec. 13. 3p . the contribution of 
infinitely large wave vectors is unimportant here. Hence, calculations here and in Sec. [3] yield 
consistent estimates of the renormalized Hamiltonian. 



5 Renormalization up to the order of e 

5.1 Vertices of the order O (e^) 

According to our diagrammatic representation, vertices of the order O (e^) in the renormalized 
Hamiltonian are those made by coupling three original vertices of ([2]) , as well as other diagrams of 
such topology. The following rigorous statements are used to find the topological pictures, made 
by coupling lines, corresponding to the vertices of the order O (e^). These are the topological 

pictures obtained when the dashed lines of the vertices <^ shrink to points, so that we 

deal with instead. 

Lemma 1. Any connected diagram made of three vertices ^^>\ contains a subset of two 
coupling lines connecting the nodes as in the diagram — | — | — | — . 

Proof Let number the vertices as vertex 1, vertex 2, and vertex 3. Since the diagram is 
connected, vertex 1 must be connected at least with one line either to vertex 2 or to vertex 
3, and these two connected vertices must be connected with at least one line to the remaining 
vertex. Hence such a subset of connecting lines as in the diagram — | — | — | — always exists. □ 

Lemma 2. Any diagram of the renormalized Hamiltonian containing a subgraph (which is 
not the whole diagram) of topology — • — , — Q— , , 3 , ^ , etc., i. e., 
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a subgraph with two outgoing lines, can be nonvanishing only if the lines of this subgraph are the 
internal ones of the whole diagram. 

Proof. The considered subgraphs have opposite wave vectors k and — k related to their two 
outgoing lines, as it follows from the fact that the sum of the wave vectors entering each node in 
a Feynman diagram is vanishing. On the other hand, by definition of the RGT, the internal and 
the external lines in the diagrams of the renormalized Hamiltonian always have different values 
of I k |. Hence, the diagram can be nonvanishing only if both outgoing lines of the subgraph are 
either internal or external lines of the whole diagram. Only the first scenario is possible, since 
the subgraph has to be connected to the rest of the diagram. □ 

Lemma 3. Any connected diagram with at least four external lines, which is a nonvanishing 
diagram of the renormalized Hamiltonian and is made of two vertices and one vertex 

— • — , contains a subset of two coupling lines connecting the nodes as in the graph — | * | — . 



Proof. The lines of the vertex — • — have to be the internal lines of the diagram according 
to Lemma 2. Hence, there are two possible ways how these lines can be coupled to the rest of the 

diagram: (1) they are connected to two lines of one of the vertices as ; and (2) 

they are connected to both vertices as — | * | — . In the first case, the remaining vertex 

X has to be connected as ^ to obtain a dia gram with at least four external lines. 
This diagram, however, is vanishing according to Lemma 2. Hence, only the second possibility 
remains, i. e., the subset of coupling lines contained in the diagram — | * | — always exists. 
□ 

It is easy to verify that nonvanishing connected diagrams made of two vertices — • — and 
one vertex can contain no more than two external lines. Thus, Lemma 1 and Lemma 3 

imply that the diagrams of order O (e^) with at least four external lines can have only the 
following topologies: — | — [— | — , — | * | — , and those ones obtained by coupling the lines 
in these graphs. Using also Lemma 2, it allows to find easily all the nonvanishing diagrams of 
this type. In such a way, the renormalized Hamiltonian has the form 

H = + ^(4) + ^(6) ^ ^(8) ^ (^^4^ ^ (72) 

where ff^^™-) includes all the (/j^^-type vertices representable by the diagrams with 2m external 
lines. The corresponding algebraic representation is 

^(2m) 

__ ^ yl-m Qm(kl,qi,...,km,qm)(/9j^_kiVii,qi •••'/'im,k„y'j„,q„ • 

(73) 

In such a form Qm (ki, qi, . . . , k^, q™) vanishes unless J2i ki + qi = 0. The dependence of these 
weight factors on various parameters is not indicated here. 

As regards the quartic part H^'^\ we need to separate the contribution corresponding the 
ordinary vertex, which is provided by the constant part of the weight function Q2 in (|73|) 
evaluated at kj = q^ = 0. For this purpose we can represent each diagram of as 









r 










^ + 





X X • (74) 
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In this symbohc notation 



is any diagram with four external hnes, whereas 
is the same diagram from which the contribution corresponding to the ordinary ip'^ vertex is 



subtracted. The latter one is represented by the last term in (j7i|) . where q I I q is the 
diagram with amputated external lines having zero wave vectors. A particular example 



>::::::< = >■■/■■■<. + o<::> x x (75) 

refers to the diagrams already considered in Sec. [3l A problem here is that not all diagrams of 
are convergent when the integration region for the internal lines [A, A^] is 



the kind 

extended to infinity at C — >• oo. To avoid possible unphysical divergence and instability of the 
RG flow, a suitable representation of the renormalized Hamiltonian should be used, where all 
vertices are convergent at (" — )■ c«. It is reached by making appropriate zero-k subtractions not 

only from the entire ip^ diagrams, but also from their internal parts, like, e. g., 

and • In the latter case, it is not necessary to make an extra subtraction of the 

zero-k contribution from the whole diagram, since such a contribution is vanishing. Appropriate 
subtractions should be used also for vertices of higher than ip'^ order to ensure their convergence 
at C — oo- 

5.2 The renormalized Hamiltonian and RG flow equations 

According to the analysis of Sec. [5l an appropriate representation of the renormalized Hamilto- 
nian is 



H 

T 



1 



^ (r + ck' + 0(k)) I P +1 >--< + a« S ( >X< 




+ a. 



i^^S > 



,(2) 



>••-•< )+a8E(4-|-f)+0(. 



(76) 



where S(-) denotes the sum of all such diagrams made of vertices 



and 



which yield the given picture when the dashed lines shrink to points. Besides, the combinatorial 
weight coefficients, including the n~dependent factors, are normalized in such a way that their 
sum is 1 at n = 1. The relations of the kind 



oO X 

oO X >^ 



(77) 
(78) 



have been applied to reduce the number of different kind of vertices in the representation of the 
Hamiltonian. The wave vectors of dotted lines are within k G [A, A^] with C being updated as 
C = s C under the RG transformation Rg . 
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The other Hamiltonian parameters are updated as follows: 



u 



n + 8 



oO + 



+ 6n + 20 



-u 



( o<:> 



+ (5n + 22)nM o< + 0<>r + 2 0< =. + 2 o\ =. + OV 



+ 



(n + 2)(n + 8) 



„<B> 



+ 0- 



,(4) 



,(3) 



,(4) 



^ u^s^e (l - (n + 8)u oO )+0(^£ 
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+ (n + 8)nV o<^ 
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32 
27 



-iti^s^' (l - (n + 8)ii oO ) + O (e") 
iuV + O (£*) 

A„3 + 0(,4). 



^ '2 I 

:W + 



0{e') (79) 

(80) 
(81) 
(82) 
(83) 
(84) 
(85) 
(86) 
(87) 



r + oO ) + (e^) is the updated 



Here k G [A/s, A] holds for solid lines, and r' = , ^ 

value of r in the lowest order of the e-expansion. We have skipped the refined RG fiow equations 
for r, c, and 9{k), since those already considered in Sec. 13. H are sufficient for our further analysis. 

These RG flow equations are proven by direct veriflcation, as in Sec. 13.21 Besides, due to 
the convergence of all vertices of (|76p at ^ — )• oo (as it is verifled by estimating the k-space 
integrals like in Sec. 13. 3|) . this renormalization is consistent with such one, where the entire 
interval k G [A, A^] is split in slices like in Sec. [H 

The fixed-point values of the parameters a"l^ , a^^ and ag are trivially related to u* and r*. 



Namely, (^a 



(1) 



-(9/16)n*2 + o(e3), (4^))* = (9/8)m* V + o (e^), (^i))* = -(l/8)n* 



+ 



o (e ), and similarly for other parameters. The value of u* is 



2c2A=e 



X (l + ai{n)e + O (e^ 



{n + %)Kd 

The coefficient ai(n), calculated from (f79l) . (f56]l . and is 



ai{n) 



4(5n + 22) 



Q 



n + 2 
(n + 8f 



(90) 
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Here Q = (5(s)/lns, where Q is given by the diagram expression 



Q{s) 




+ 0< 



+ 2 0' 



+ 2 o-:. 



+ 



(91) 



evaluated at d = 4 and c = 1. The quantity Q{s) is proportional to Ins, since Q{siS2) = 
Q{si) + Q{s2) holds at ^ — )• oo as a diagram identity, which can be verified by the method 
outlined in Sec. 13. 4[ Hence, Q is s-independent and, therefore, u* is s-independent at least up 
to the order, as expected. It is convenient to consider the limit s — )• oo to determine Q from 



Q 



lim 



lim 



lim < 

s— >-oo 



Q{s) 

In s 

1 
Ins 



1 
Ins 







= lim < 






In s 


...0 






o\ ; ; - 





.0 



Ki 




1 

+ 2 



Kl 



(92) 



The wave vectors are rescaled to € [1, s] (i. e., this interval corresponds to the dotted lines in 
the diagrams of (j92p ) and the zero-k contribution of the internal diagram block is separated in 

the second line of The large-A; asymptotic estimate for k*-/i^-' ) obtained in Sec. 13.31 is 
used in the third line of ([92]) . Inserting this result into ([90]) . we obtain 



ai{n) 



3(3n + 14) 
(n + 8)2 



(93) 



Summarising the results of this section, one has to note that, at the order, the renormalized 
Hamiltonian contains several extra vertices as compared to the initial or bare Hamiltonian. All 
these terms, including the ip^ vertices and the ip^ vertex, are relevant in the sense that they do 
not vanish at the fixed point. 



6 Expansion at a fixed spatial dimension d 

Apart from the e-expansion, we have tested also an alternative approach, where the coupling 
constant u is considered as an expansion parameter at a fixed spatial dimensionality, i. e., fixed e. 
In this sense the discussed here approach is similar in spirit to that widely used for calculations 
in three dimensions, d = 3, as described in [Ml [HI [TU [iTl [18]. This known approach uses the 
Callan-Symanzik equation instead of the Wilson's equation ([3]). However, if both are correct RG 
equations, then they should provide consistent results. Here we use ([3|) to make the expansion 
in powers of n at a fixed but small e. 

In this method, all the rescaling factors of the kind s^ have to be retained in their original 
form without the expansion in powers of e. Following the approach of |14[ [T5l [T6l [T7[ llSj . 
the critical exponents should be expanded in powers of the coupling constant, estimating their 
universal values at u = u* . Therefore, the correction factors like s~'^ can be omitted in the first 
approximation, since one finds that rj = O (u^)- Using this idea, and following the calculations 
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in Sec. 13.11 we determine the fixed-point values of u and r in the lowest order of the theory, 
i. e., u* and f*, as 

KdS^ n + 8 ' ^ ^ 

r = — TT, e • (95) 



These values are obtained as the fixed-point solutions of (I2ip and (1221) . retaining the relevant 
terms without their expansion in powers of e. In this case the terms up to the lowest-order 
diagrams, o"^^^ iii (EH) and qO in ([2^ . have been included. 

As we see, the fixed-point values depend on the scale parameter s, thus pointing to an 
internal inconsistency of this method. It is in contrast to the results of the e-expansion, where 
such a problem was not detected. 



7 The two— point correlation function 

7.1 The diagrammatic representation and scaling 

Consider now the Fourier-transformed two-point correlation function G(k) = (| y?j^k P)- Its 
diagram expansion contains all connected diagrams with two external lines, having wave vectors 
k and — k. Thus we have 

G(k) = ^ + + ^-^^ +... = ^ , (96) 

1 — 2Go(kj2j(kj 

where — - represents the sum of diagrams of this kind, which are made of vertices of 
—H/T and are irreducible in the sense that they are not representable as a chain of such blocks. 

Here S(k) = ^^-X- ^ X ^ is such irreducible diagram block with amputated external lines. The 
quantity — 2E(k) is known as self-energy |2j. The well known expansion (|96p leads to the 
equation 

X(k) = ck"^ - 2S(k) , (97) 

where X(k) = 1/G(k) is the inverse of the two-point correlation function, and G'o(k) is set 
to l/{ck'^) within the e-expansion. In this case, the diagrams contain also insertions of the 
second-order vertex {l/2)J2{f + ^(k)) | V3j,k P- The wave vectors within q G [0, A] are related 

to the solid coupling lines between different vertices in these diagrams. The vertices themselves 
can contain internal dotted lines with q £ [A, A^], produced by the RG transformation, ( being 
the renormalization scale in the RG flow equations. 

As it is well known, the two-point correlation function rescales in accordance with 

X{k; RsH) = s^^'''X{]i/s; H) (98) 

after the RG transformation Rg , transforming the original Hamiltonian H into the renormalized 
one RgH. It is a general exact relation, trivially following from the fact that the integration 
over the Fourier modes with A/s < /c < A in the flrst step of the RG transformation does not 
alter the correlation function for k < A/s, whereas the s-dependent factors in (I98p compensate 
its rescaling in the second step of the RG transformation (see Sec. [2|). 
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7.2 The two— point correlation function at the fixed point 

In the following, we will examine the two-point correlation function (or its inverse) at the fixed 
point up to the order. According to (198p . X(k) = ak'^~^ holds within k < A for the fixed-point 
Hamiltonian H* (since RsH* = H*) with some k-independent constant a, which can depend 
on n, e, and Hamiltonian parameters. It is consistent with (I97p if 

ak\l - ^Infc) = ck^ - ^n*2 | k(|) + k0 + 3 ( k0 + k0 ) } + O {e^) (99) 

holds for A; < A, since S(0) vanishes at the fixed point. Note that q S [0, A] refers to the solid 
lines and q G [A, A^] at C — ^ oo - to the dotted lines in ([U^ . The diagrams can be evaluated at 
d = 4 in this case, and we can set c = 1, as it gives merely a common factor. By definition of 

kV.'l? , we have 

k(:|> = k0 - k'RiC) , (100) 

where 

f 1 ••■/• 1 

R{C) = lirn^ I ]^::p \ = — ^ In C + const at d = 4, C ^ oo , (101) 

as consistent with the already considered estimation of such term (after rescaling the wave 
vectors from q G [A, A^] to g G [1,C]) at the end of Sec. 13.51 Inserting (jlOOp into (f99|) . and 
summing up the diagrams, we obtain (at c = 1) 



ar(l -r/ln/c) = r —u*^f{k,C) + 0[e'^) at C ^ oo , (102) 

where 

/(^,C)= k0 -fc'ii(C) (103) 

with g G [0, A^] corresponding to the dashed lines. Using (jlOip and rescaling the wave vectors, 
we find that 

/(A:)= hm f{k,0 (104) 

C— >-oo 

obeys the equation 

f{sk) = (^f{k) + ^Kj k^ In (105) 

with the solution 

f{k) = Bk^ + ^Kl k'^lnk , (106) 

where i3 is a constant. According to this, Eq. ()102p is indeed satisfied with appropriate a = 
a{n,e) at r] given by ()64p . In such a way, we have verified that the e-expansion of X(k) at the 
fixed point is consistent with the scaling Ar(k) = ak'^~^ within A: < A up to the order of e^, at 
least. Besides, such analysis provides one more method to determine the critical exponent r]. 
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7.3 The two— point correlation function on the critical surface 
7.3.1 Some consistency relations 

The critical correlation function diverges at /c — >• 0, so that ^(0) = and S(0) = must hold 
in this case, which means that (j97p becomes 

X{k) = ck"^ - 2[S(k) - S(0)] . (107) 

It is easy to verify that the condition S(0) = is indeed satisfied in the lowest order of the 
e-expansion, where the equation of the critical surface reads 

r = rc{u) = oO +0 (e^) (108) 

with k G [0, A] for the solid lines. The condition S(0) = should hold on the critical surface 
also at higher orders of the e-expansion. We have tested this expected property up to the 
order and have not revealed any contradiction. 

According to (jl07p . the equation for X(k) on the critical surface reads 

^(k) = ck' - ^^u' I k(|;) + k0 + 3 ( k0 + k0 ) } + O {e') . (109) 

Based on this equation, we will consider the rescaling of X(h) in the renormalization process, 
starting with the bare Hamiltonian at {u — u*)/u* = do and c = cq and performing certain 
number ni of RG transformation Rg, in such a way that (n — u*)/u* becomes equal to 6, where 
6o and 6 are small constants. In this case, the fixed-point value u* of the coupling constant is 
determined at the current value of the parameter c, so that it is a slightly varying quantity due 
to the renormalization of c. Such a choice is convenient, since we do not need to determine the 
fixed-point value of c after infinite number of RG transformations. In the actually considered 
truncation of the e-expansion, the renormalization of u is given by ()2ip . from which we evaluate 
m = m{6o, 5, s, e). In this case, m ~ 1/e is determined with uncertainty of 0(1), so that it can be 
rounded to appropriate integer value at any positive s — 1 of order unity. The renormalization 
of {u — u*)/u* can be performed at integer m for e-independent <5o and 6 up to any higher 
order of the e-expansion by adjusting s as s = sq + sie + S2e^ + • • • . The entire renormalized 
Hamiltonian after m ~ 1/e RGT is well defined within the e-expansion, since we always can 
reach the required accuracy, truncating the RG fiow equations at high enough order. 

Our aim is to compare the perturbative rescaling of X(k) with the exact one ()98p at /c — t- 0, 
taking into account the expected asymptotic form of Ar(k) on the critical surface. Namely, for 
the initial Hamiltonian we have 

X{k;Ho) = A{6o,s)k^-'^ i^l + Y,ai{6o,e)k'^' +J2ai{6o,e)k'^''^ at k^O, (110) 

where the argument Hq means the initial Hamiltonian, Ui are those correction-to-scaling expo- 
nents, which tend to zero at e — )• 0, whereas cDj are those ones, which tend to positive constants 
in this limit. Eq. (|98p then yields 

X(k; Rs^Ho) = A{6o,e) k^'i 1 1 + II «*('^o, e)^-™^' A:"^' + ^ a^{5Q,e)s~''^^^k^^ | at ^ , 

(111) 
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where Rs^^Hq denotes the renormahzed Hamiltonian after m RG transformations Rg. The well 
known result of the e-expansion states that the leading correction-to-scaling exponent is 

oj — oji = e + o(^e'^^ . (112) 

In the following, we will formulate as a theorem an important result, stating certain consis- 
tency relations for the coefficients and exponents in (jllOp . 

Theorem. // the perturbative e-expansion-hased renormalization (repeating the RG 
transformation Rg) o/X(k) = 1/G(k) on the critical surface is consistent with M1U\) and Mll\) 
at any small e-independent 60 and 5, such that < 6/60 < 1, defining the initial (Sq) and the 
final (5) value of {u — u*)/u* (with u* determined at the current c), then 

ooj = 3e + 0(e^) : j>l, (113) 
A{5o,e) = co + 0(e), (114) 

a,(<5o,£) = hje[-^^\o(e^) : j>l (115) 

hold, where cq is the initial value of c, and the coefficients bj are found from the equation 

6 



r] 



+ ^"T.jbA-^A =-v{l + Sf + 0{e') , (116) 



requiring the consistency of both sides at each power of 5 in the small-6 -expansion. 

Proof. Let us first determine the number m = m{6o, 6, s,e) of RGT required for transforma- 
tion of {u — u*)/u* from 60 to 5. From (|2T]) we obtain the updating rule for 5i = {ug — uD/u^, 



where ui is the value of u after ^-th RGT, and u| = n*(Q) is the fixed-point value of u determined 
at the current c after ^-th RGT, i. e., at c = q. Here we take into account that q+i = 
C£ + (e^) holds according to (p3]) and therefore u}_^_i = Uf, + (e^) fohows from ([89]) . Denoting 
xi = In \ 6e \ , Eq. (|117p reduces to 



xe+i - Xi 



-e (l + sign(5)e^^)lns + 0(e2) . (118) 



The variation of X£ becomes quasi-continuous in the limit e — )• 0, so that xi can be considered 
as a continuous function x{i), and (jllSp transforms into the differential equation 

"^""^^^ - -e (1 + sign(,5) e^) Ins + O (e^) . (119) 



d£ 

Integrating this equation, we find the necessary number of RGT 

rlT 1 r A„ /' 1 -I- /I 1 

+ (120) 



1 /■ dx , 1 , 

m = / : — — — + 0(1) = In 

ems J 1 + sign(oj ems 



5o(l + 5) 



where xq = In | 5o I is the initial and = In | 5 | is the final value of In | (5^ | . It yields 

= i±ia X + 0{e) . (121) 

Oq 1 + 
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Inserting ()12ip into we obtain 



X{k;R,..Ho) = A{do,e) fe^"" <^ 1 + ^ai(,5o,£) 



LOi/e 



(122) 

where the terms with Wj are absorbed in the last remainder term, where A > and w is the 
smahest exponent among Wj. This term is irrelevant in the e-expansion. 
From (jl07p we obtain 

X (k; Rsm Ho) = ck'^ - "^^^ [u* {1 + 6)f (^Bk'^ + ^Ki k^ In k^+o(^e^) . (123) 

The diagrams of (jl07p are evaluated in (jl23p as in Sec. 17.21 (cf. Eq. (jl06p ). since the difference 
between the cases C = oo and ( = s"^i^o,s,s,e) ^^j^^j-^ q £ \Q f^j- dotted lines) is irrelevant 
within the e-expansion. Since c is renormalized by O (e^) in one RGT according to (j23p . the 
renormalized c value after m ~ 1/e RG transformations in (I123p is c = cq + 0{e). 

The consistency with (jl22p (where ^ at e — 0) implies that X(k;RsmHQ) can be 
expanded as 

Xik; Rs^Ho) = k^ (Bo{5o,5, e) + Bi{5o, 5, e)\nk + B2{So, 6, e)(ln kf + -- •) (124) 

at e — 7- 0. On the other hand, the diagrammatic perturbative equation for X(k; Rs^Hq) on the 
critical surface can be represented in general (up to any order of the e~expansion) as 

X(k;R,r^Ho) = cF{k,6,e) , (125) 

where c = c{6Q,6,e), whereas F{k,6,e) contains only integer powers of 6 in the small-(5- 
expansion. It is because the critical value of r, as well as all Hamiltonian parameters can 
be expanded in integer powers of n = u*{l + 6), where u* = n*(c) = cPu*{l) holds with Sq- 
and 5-independent u*{l), and the diagrams associated with contain 2£ — 1 internal coupling 
lines giving the factor c^~^^. These are consequences of the general structure of the RG flow 
equations. Hence, the small-J-expansion of the ratio Bi{5o,5,£)/Bo{6o,6,s) in (I124p contains 
only integer powers of 5 for any i > 1. Assuming that coj = dj e + O (e^) holds with noninteger 
positive dj for some j, we arrive at a contradiction, since then the e-expansion in (I122p pro- 
duces Bi{5Q,5,e) /Bo{5Q,5,e) containing noninteger powers of 5. The contradiction is obtained 
also if we assume the existence of ujj of a higher order than ©(e), since the e-expansion of 
i3j((5o, (5, e)/i3o((5o, (5, e) in ()122p produces powers of In 5 in this case. We can consider also a 
possibility that ojj ~ e^^ with noninteger fij holds for some indices j. It does not lead to the 
consistency, since the expansion of i^(k, 5, e) contains only integer powers of e, coming from the 
e-expansion of u* and diagrams evaluated at c = 1. In such away, (I122p can be consistent with 
the diagrammatic perturbative equation only if (jll2p and ()113p hold. 

X{k] Rsr^Ho) = ck'^ + O (e^) = cqA;^ + 0{e) holds according to (fT23]l . Besides, the coefficient 
at In/c in ()123p is of order O (e^), and it is independent of 5q at this order. Eqs. ()114|) and ()115p 
must hold to satisfy these relations in (|122p . performing the e-expansion. Finally, (jll6p states 
the necessary consistency relation for the coefficient at e^lnA; in (I122p and (I123p . taking into 
account that [{n + 2) / {Sc^)]Klu*'^ = cr] + (e^) holds according to §5^) and □ 
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7.3.2 A further test of consistency 



Based on the relations stated in the above theorem, here we perform some test for single RG 
transformation Rg on the critical surface, starting with small {u — u*)/u* = 6q, defined as before. 
We consider the asymptotic small-k estimate of the inverse correlation function defined as 



(126) 



obtained by neglecting the irrelevant (at /c — )• 0) terms with exponents cjj in (jllOp . Follow- 
ing (j98|) . it rescales exactly as 



(127) 



after the RG transformation Rs. According to (jllSp . we can write 

5n 



ai((5o,e) = bie 



i > 1 



(128) 



where £i{5o) are some unknown coefficients. Further on, we consider the expansion up to the 
order for the ratio X"''^(k;RsHQ)/X"''^{]<.;H()), in which case these unknown terms cancel, i. e.. 



:i 



5.? 



+ e 



A°^(k;fi,go) 
A-(k;Fo) 

(129) 

The sum over i in (|129p is meaningful at least as the small-[5o/(l + 5o)] expansion, which gives 
the stated here final result after re-expanding in powers of 6o in accordance with ()116p . Thus, 
this is the result which holds in all orders of the small-Jo expansion. 

In the following, we calculate the ratio X^^{\<.; RsHQ)/X'^^{]i.; Hq) from the diagram equa- 
tion ()107p and compare the results. In this case we need to include those terms, which correspond 
to A°^(k; i?si7o) and X"'^(k;HQ) in the e-expansion. Thus, the terms containing k'^{lnkY with 
£ > must be included, whereas those containing A:'^(lnA;)^ with > 2 have to be neglected. 
Any quantity endowed with superscript "as" is calculated in this way. Furthermore, both 
A"*(k; RgHo) and A"**(k; Hq) are proportional to the initial c value cq, so that we can calculate 
their ratio at cq = 1 without loss of generality. By definition, the term ^(k), represented by 

kV'l'.' ) tends to zero faster than ~ A;^ at — >• 0, and therefore is neglected. Thus, the diagram 
equation at cq = 1 yields 



A-(k;i/o) 
X'^'ik-RsHo) 



k —u [l + do 



ck 
3 



(l + (5o) 







+ 



+ 

as 

+ 



(130) 



(131) 



where q E [0, A] holds for the solid lines and q S [A, As] - for the dotted lines. Here we have 
set u = n*(l + 6q) in both expressions, since u'^ after the RGT is varied only by O (e^). The 
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quantity ^ ^ cancels in X^^(k;RsHo)/X"-^(k;Ho), calculated up to the order. The 

asymptotic estimate for the sum of two other diagrams at d = 4 and c = 1 is 



+ 



A;2$(0,s) 



where 



^(Ai, A2) := lim 



1 



A:-s>0 1 fc2 



(132) 



(133) 



with q £ [Ai, 1] for the solid lines and q G [1, A2] for the dotted lines. Eq. (jl32p is obtained by 
normalizing the wave vectors {q — )• q/A) and extracting the relevant (at — )• 0) contribution 
oc /c^, neglecting corrections like ~ /c^, ~ k^, etc. The renormalized c in p3ip . calculated 
from (f23]l . is 



1 



, _ !i±l„.=(i +,,,2 {i. ( + O (134) 



with q € [A/s,A] corresponding to the solid lines in the diagram evaluated at c = 1. If this 



diagram is completed by 3 k' 



.0 



+ kt 



with q E [A/s, A] for the solid lines and q G [A, cxd] 

for the dotted lines, then the obtained expression lim i -r^- (•) > is exactly oc Ins, as it follows 

from the diagram identity (j62p . The proportionality coefficient is —K^/A, in accordance with 
the calculations at s — )• 00 in Sec. 13.51 Hence, rescaling the wave vectors, we can write 



lim , , „ 



Ins - 3$(l/s,oo) 



(135) 



d=4 



with <I>(l/s,oo) defined by (jl33p . Summarizing all the derived here diagrammatic relations for 
X'^''(k; RsHq) and X"*(k; Hq), and also [{n + 2)/(8c3)]K|n*2 = cr] + O (e^), we finally obtain 



X"*(k; RsHq 



1 + V 



(1 + So)^ - 1 



In s 
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It is evident that (I129p and (I136p agree with each other if and only if 

(j)(s) = <^{0,s) - $(l/s,oo) = 



(136) 
(137) 



holds for s > 1. 

Consider now the vicinity of s = 2. It is a special point for one diagram of (/){s), i. e., k{^ 
with q E [0, 1] for the solid lines and q G [1, s] for the dotted line. Some singularity is expected 
here, because the integration over the region with q > 2 for the dotted line gives vanishing result 
at A; —7- and, therefore, the contribution of the mentioned diagram becomes s-independent at 
s = 2, when s is increased. By definition, this contribution is equal to 0(1, s), where 



(/.(Ai,A2 



F-5o i 



d=4 



(138) 



with q S [0, 1] corresponding to the solid lines and q G [Ai, A2] - to the dotted line. Using the 
diagram identity 

0(1,00) = 0(1, s) + 0(s, 00) , (139) 
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we will evaluate the singular part of (/>(!, s) by calculating 

0(s,oo) = lim|l^ |[/(|q+k|)-/(g)]g-Vg| , (140) 



where 



/(g) = qO =-i-4 1 A;-2|q + kr2e(l-|q + k|2)d4A: 



k>s 



oo 



2K4 / , „ /■ G (1 - g2 _ 2g/c cos e-k"^) 2.,. . ^ 



TT J J g2 + 2qk cos 9 + k'^ 

s 

Here A; G [0, 1] corresponds to the solid line and k £ [s, 00] - to the dotted line in the diagram. 
For 1 < s < 2, we obtain /(g) = if g < s — 1, and 

2/r4 /• , „ f sin^ede 



/(g) = — - kdk / — -r (142) 

-K J J g2 + 2g/c cos 6* + A;2 ^ ' 

with ^(g, k) = arccos [(1 - g^ - A;2)/(2gA;)] if s - 1 < g < 2. 
In the following we use the small-k expansion 

/(I q+k|) =/(g) + — Ucos^ + ^sin^ej +-^A;2cos20 + o(A;3) , (143) 

where 6 is the angle made by q and k. Inserting this into (I140p . we obtain 

<^(s,oo) = ^ (2/(1) + {di/dq)\^^^ (144) 

for 1 < s < 2, taking into account that /(g) and dl{q)/dq vanish at g < s — 1. Using (11420 . we 
find that the leading singularity 

(^(s,oo) ~ ^ A^/^ A = 2-s^0 (145) 

is provided by the derivative term. According to (|139|) . —</>(!, s) has exactly such singular part. 
Furthermore, since s = 2 is not any special or singular point for the diagrams of (p{s) in (jl37p . 
except only the one related to 4'{l,s), we conclude that (/)(s) has the same singularity as (p{l,s). 
Hence, 4>{s) = cannot hold, and the consistency condition (jl37p is not satisfied. 

This fundamental contradiction implies that the perturbative e-expansion-based renormal- 
ization of X{h) on the critical surface is not fully consistent with the one produced by the exact 
rescaling (j98p of the expected asymptotic expansion (|110p at small k. Note that the renormal- 
ization considered in Sec. 17.3.11 is an infinite-scale (m — )• 00) renormalization at e — )• 0. If the 
e-expansion works correctly there, then it turns out that an inconsistency appears on a finite 
renormalization scale considered in this section. The obtained here contradiction shows that the 
e-expansion fails to give correct results in the actual test, despite of the fact that it is formally 
well defined and stable. One has to note that the contradiction, apparently, is not caused by 
the violation of the rescaling rule (j98p in the e-expansion, but rather by the inconsistency with 
the asymptotic expansion (jllOp . rescaled according to (I98p . Indeed, we have not obtained any 
contradiction by testing the consistency with (j98p alone. 
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7.4 The two— point correlation function above the critical point 

Consider now the inverse of the two-point correlation function above the critical point. It implies 
that X(k) has to be positively defined for all k when the diagrams are represented by k-space 
integrals. A problem here is such that the equation contains divergent diagrams. This difficulty 
can be overcome in a modified model where the lower cutoff A of the wave vectors is introduced. 
In this case q G [A, A] corresponds to the solid coupling lines in Eq. (|97p . and the e-expansion 
is well defined for < A < A. Namely, we have 

X{k) = ck' - ^^n' I k<|) + k0 + 3 ( kG + k(/) ) } 

+ (r + ^noO)(l-^no<:>)+o(e^) (146) 

with q E [A, A(^] for the dotted lines at the given renormalization scale C- Here C = 1 corresponds 
to the bare Hamiltonian. Note, however, that A is not constant, but is updated in the RG 
transformation Rg as A' = sA, so that A = AqC holds, where Aq is the initial value of A. Hence, 
it is meaningful to consider the renormalization scales 1 < C < A/Aq. We have verified that the 
rescaling property ()98p is satisfied up to the order in this model with r = 0{e) and u = 0{e). 
On the other hand, we obtain certain paradox: the diagram o^^^^ diverges and, therefore, 
the e-expansion becomes unstable at Aq — )• despite of the fact that the correlation function 
in the high-temperature phase should not be singular even at Aq = 0. A standard idea, widely 
applied in the perturbative renormalization, is to subtract divergent terms. One believes that 
a physically correct result is obtained by such a method, also called renormalization. In the 
modified model, however, all terms are finite, i. e., there is no reason to subtract any of them at 
finite values of Aq. Another question is whether a correct result is obtained in the original model 
with Aq = 0, if the divergent diagram o^^^ is simply discarded. The answer is negative. 
Indeed, it can be easily checked that ^(k) obtained in this way is inconsistent with the exact 
rescaling relation ()98p . The latter finding is not surprising: in the diagrammatic perturbation 
theory this relation, similarly as the semigroup property and the s-independence of the fixed 
point of the RG transformation Rg, is a consequence of certain diagram rescaling and summation 
rules. These rules cannot work to give the expected result if some diagram is simply omitted. 

8 Discussion 

There are a lot of perturbative RG studies of the 99^ model made in the past (see [21 [U HI [19] for 
a review). They can be classified as approximative treatments of the perturbation theory in view 
of our analysis, since a set of apparently irrelevant terms is neglected. Although a perturbative 
RG approach to critical phenomena never could be considered as a rigorous method, it is, 
nevertheless, possible to calculate exactly the perturbation terms and to take them into account 
in a systematic way. In particular, we find that there are many different terms to be included 
in the renormalized Hamiltonian up to the order of e^, and these terms are relevant in the sense 
that they do not vanish when repeating the RG transformation unlimitedly many (m — ?■ 00) 
times. To the contrary, one usually finds that only few terms, i. e., those included already in 
the initial Hamiltonian ([2]) are relevant. A confusion about this might be caused by the fact 
that in the usual (nonrigorous) treatments one looks for the behaviour of individual terms and 
finds that most of them are shrinking in the renormalization procedure. It, however, is not a 
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rigorous and even not a valid argument, since the ff^^*") parts of ()72[) (as well as any part of 
^(2m) represented as a sum over diagrams of certain topology) result from many such individual 
terms, summing up to yield relevant contributions. 

The results of perturbative calculations up to the order of are reported in literature. These, 
however, are based on the Callan-Symanzik equation, but not on the Wilson's equation ([3]). 
The Callan-Symanzik equation, in fact, represents a scaling property for the Hamiltonian of the 
form ([2]). The approach, thus, relies on the assumption that all the relevant terms are included 
in the initial Hamiltonian, only the existing here coupling constants being renormalized. This 
assumption is not supported by our analysis: the perturbative renormalization based on first 
principles (i. e., avoiding such assumptions via using ([3])) produces additional relevant terms. 

One believes that the e-expansion, combined with the perturbative RG transformation, is 
able to describe correctly the small-k asymptotic of the correlation functions and the related 
critical phenomena. The contradiction derived in Sec. 17.3.2] causes serious doubts about it. Due 
to the formal character of the e-expansion, one might happen that it describes the behaviour 
within any fixed range of wave vectors at e — )• rather than the fc — )• limit at a small, but 
fixed e. Indeed, the contradiction might be caused by the fact that the non-Gaussian part of 
the Hamiltonian cannot be considered as a small perturbation, like in the e-expansion, if one 
considers just the small-k contribution responsible for the critical fluctuations, since the Fourier 
amplitudes of the relevant ¥'j(x) configurations diverge when A; — )• at criticality. 

The problems and paradox discussed in Sec. l7.4l rise a question concerning that how the two- 
point correlation function slightly above the critical point can be correctly calculated within 
the e-expansion, and whether a correct calculation method does exist here. In particular, it 
turns out that a simple discarding of divergent diagrams is not a correct approach, since it 
leads to the violation of the rescaling rule ()98p . In fact, certain method of calculation of the 
two-point correlation function has been considered in [12], where the discussed here paradox 
does not appear. It, however, is not based on the e-expansion, but on certain grouping of 
diagrams, in such a way that the true correlation function instead of the Gaussian one is related 
to the coupling lines. It naturally ensures the convergence of the diagrams of the reorganized 
perturbation theory and allows an analytic continuation of the obtained self consistent equations 
from the region r > to the vicinity of the critical point. 

9 Conclusions 

1. A diagrammatic formulation of the perturbative renormalization has been provided (Sec. [2]), 
which makes calculations of perturbation terms straightforward and transparent, avoiding 
any intermediate approximations. 

2. The RG flow equations, including all terms up to the order of e^, have been considered in 
Sees. [3] and m The tests of the expected properties, such as the semigroup property and 
the existence of an independent of the scale parameter s fixed point, have been performed. 
These properties are satisfied up to this order of the e-expansion. No inconsistency has 
been found at the e^ order, as well (Sec. [5]). However, some internal inconsistency has 
been detected in an alternative approach, where e is fixed and the coupling constant u is 
an expansion parameter (Sec. [6]). 

3. At the e^ order, the renormalized Hamiltonian contains several extra ((/9^, ^p^ , tp^) vertices, 
as compared to the bare Hamiltonian. All these terms are relevant in the sense that they 
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do not vanish at the fixed point. 

4. The two-point correlation function has been considered in Sec. [7] within the e-expansion 
up to the order of e^. It is consistent with the expected power-like behaviour at the fixed 
point (Sec, \7.2\i . Nevertheless, a fundamental contradiction with the exact rescaling of the 
small-k asymptotic expansion on the critical surface exists (Sees. 17^ I7.3.2p . 

5. Problems and certain paradox, related to the e-expansion of the inverse correlation func- 
tion slightly above the critical point, have been considered in Sec. 17.41 

Concerning the last three points, several related issues have been discussed in Sec. [HI 
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